A collective modulation instability of multiple four-wave mixing 
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We investigate the modulation instability of multiple four-wave mixing arising from a dual-frequency pump 
in a single-mode fiber or waveguide. By applying the Floquct theory on account of the periodic nature of 
four-wave mixing, we reveal a collective type of instability occurring in the anomalous dispersion regime. Our 
interpretation of the linear stability analysis is validated by the numerical solution of the nonlinear Schrodinger 
equation. © 2013 Optical Society of America 
OCIS codes: 190.4380, 190.4370. 



Two co-polarized pump frequencies ± (detun- 
ing 2f2fi) propagating in a Kerr medium produce multi- 
ple four- wave mixing (FWM), i.e. generation of sideband 
pairs at ujQztnfld with n odd integer. FWM is well known 
in fibers [IHl] , and is being investigated also in other set- 
tings, e.g. semiconductor photonic crystal guides [5]. De- 
spite long-lasting studies of FWM, a fundamental aspect 
such as the stability of the process has been addressed 
only recently [5], though its impact on applications (e.g., 
regeneration [7], parametric amplification [5] and pulse- 
train generation [5]) can be quite important. 

In this Letter we investigate modulational instabil- 
ity (MI) of FWM. A well-known approach to MI of 
two pumps is based on incoherently coupled nonlinear 
Schrodinger (IC-NLS) equations, entailing coupling only 
through cross-phase modulation (XPM) [TU1[TT]. While 
such model is valid for XPM-induced MI in appropriate 
polarization configurations |12I13| , in the scalar case con- 
sidered here (i.e., detuned pumps with parallel polariza- 
tion), such model does not account for FWM. Therefore 
the IC-NLS model has been argued to be inconsistent be- 
cause the predicted MI is resonant with FWM products 
n = 3 (wo ± Srid), neglected in the model [M]. 

Here we show that the MI analysis can be carried out 
by accounting for FWM, if one adopts the technique orig- 
inally developed in Refs. [T5l[T6]. Using a more general 
model, namely a single NLS equation, which is known 
to account for multiple FWM [3116], we describe FWM 
by use of four-mode truncation, assuming that the dy- 
namics of higher-order FWM sidebands is locked to the 
first four lines (the input dual- frequency pump and the 
first-order generated FWM sidebands). In this way, the 
recurrent (periodic) nature of FWM evolution [SjH] can 
be exploited to assess the stability of the process against 
the growth of additional sideband pairs around the four 
waves by calculating the Floquet exponents of the lin- 
earized problem. We interpret the gain curves calculated 
in this way and validate our interpretation by direct nu- 
merical integration of the NLS equation, finding a col- 
lective MI of FWM in the anomalous GVD regime. 

Let us consider the dimensionless NLS equation for 
a slowly-varying envelope u{z,t) of a guided mode at 
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where /3 — k" /\k"\ is the sign of group- velocity disper- 
sion (GVD), z = Z/Zni is the distance Z in units of the 
nonhnear length Z„i = (7Pt)~\ and t ^ {T - k'Z)/Ts 
stands for the retarded time in units of Tg — \J Zni \k"\, 
Pt is the total input power injected in the mode {u 
is normalized so that it carries unit power), and 7 = 
kon2i /Aeff is the standard nonlinear coefficient. 

Following Ref . [3] , we describe the dynamics of FWM 
by means of four-mode truncation: 



(2) 



u{z,t) = -^[Api(z) cxp(iilt) + Ap2{z) exp{—int) 

+^,,1(2:) cxp{i3nt) + As2{z) cxp{~i3nt)], 

where Api^p2 and Asi^s2 are the cw complex amplitudes 
of n = 1 (pumps) and n = 3 (generated sidebands) FWM 
orders, respectively, 20 = 2ndTs being the normalized 
frequency detuning (henceforth all frequencies denote ac- 
tually normalized angular frequency offsets from ujq). By 
substituting Eq. ([2]) in Eq. ([T]) and collecting terms at 
frequencies ifi, ±3Q, we obtain a set of coupled-mode 
equations which govern the evolution of Api^si- 
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whereas Ap2.s2 obey Eqs. ([3]) with subscripts 1 and 2 
interchanged. The periodic solutions with period Zp of 
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Eqs. ([3]) provide a good approximation to multiple FWM 
as long as O > 1 [3-- In order to investigate the MI of 
the periodic FWM process, we substitute in Eq. ([T]), a 
time-perturbed version of Eq. ([2]) , namely 



Akiz,t)= [Vw) +efe(^,Oj exp(i0fc(z)), (4) 

where k = pl,p2,sl,s2, rjk and 0^ are, respectively, 
power and phase of the periodic solutions Ak = 
rjf. exp(i(j)k) of Eqs. ([3|), while we assume harmonic per- 
turbations ek{z,t) = ekaexp{—iSujt)+eksCxp{iSujt), cor- 
responding to Stokes-antiStokcs pairs at frequency offset 
duj from each FWM component. Retaining only linear 
terms in Ska.ks, we obtain the following system, which 
describes the linearized evolution of the perturbation in 
a neighborhood of the dynamics ruled by Eqs. ([3]) 
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Mp (Su) X, 
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where X(z) = (Cpis 

is the perturbation column array, and Mp[dui) is a com- 
plex 8x8 matrix, which depends on 6u!, and whose 
elements are expressed in terms of periodic functions rjk 
and (/)fc. Since the matrix Mp((5a;) is periodic with period 
Zp, Eq. ([5]) lends itself to a Floquet analysis in order 
to assess its stability. Floquet theory prescribes that 
the monodromy matrix of the system is constructed by 
integrating Eq. ([5]) over the period with a set of linearly 
independent initial values (e.g. the identity matrix). 
The eigenvalues A of such matrix are called Floquet 
multipliers: the perturbation is amplified exponentially 
with gain g = 21og(|A|)/2:p, whenever |A| > 1. We 
proceed to calculate the Floquet multipliers by jointly 
integrating over z = [0, Zp] Eqs. ([3]) and the linear 
system in Eq. ([5]). Here we restrict ourselves to discuss 
symmetric FWM from initial conditions Api^p2{0) = 1, 
^,si,,s2(0) = 0), reporting detailed results obtained for 
51 = 1.5 as an illustrative example of the regime where 
FWM is non-negligible, while preserving the validity of 
the truncation {il > 1). We also assume, for the time 
being, anomalous GVD. 

The gain curve reported in Fig. [Ija) shows the exis- 
tence of three unstable branches labeled 1, 2, 3. Impor- 
tantly, the first two branches are reminiscent of those due 
to XPM-MI (dashed curves) obtained from the IC-NLS 
model [lOj . Let us remind that, in particular, band 2 ap- 
pears to be nearly resonant (Slu ~ 2il) with n = 3 FWM 
products and, as such, it was argued to be inconsistent 
with the IC-NLS model, where FWM is neglected [Ti] . 
Our result implies that such simple inconsistency argu- 
ment cannot be invoked since the band survives when 
FWM is accounted for. However, by examining the com- 
position of the unstable modes (i.e., the eigenvectors) 
in each of the three branches, we reach the conclusion 
that, while branch 1 is effective for MI, branches 2 and 
3 are fictitious in the sense that they do not give rise to 
different observable MI sidebands. The reason is as fol- 
lows. The unstable mode associated with branch 1, dis- 
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Fig. 1. (a) Gain curves from Floquet analysis (solid line) 
and IC-NLS (dashed line) vs. MI frequency Soj, for pump 
detuning ft = 1.5 and anomalous GVD {/3 = —1). (b-c- 
d) Structure of the unstable eigenvector at peak of bands 
1 (b), 2 (c), and 3 (d). The dashed red lines locate the 
FWM modes. The arrows connect the dominant MI side- 
bands to the Ak component they pertain to. 
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Fig. 2. (a,c,d) Output spectra from Eq. ([T]): anomalous 
GVD, fl = 1.5 (a) and fl = 4 (c). In (a) the gain band 1 
around pumps are reported in thin solid red for compar- 
ison; normal GVD, = 1.5 (d); (b) Periodic evolution 
{zp = 0.71) of pump power fraction relative to (a). 



played in Fig.[IJb), implies the growth of nearly symmet- 
ric components around both pump waves, and weaker 
(nearly symmetric as well) components around FWM 
sidebands. Viceversa the growing modes associated with 
band 2 and 3 turn out to have one pair of strongly promi- 
nent components (out of the four possible pairs) , namely 
£piaj £p2s (frequencies ~ Suj and — + 5uj) for band 
2, and £sia, £s2s (frequencies 3f2 — Slo and — 3fi + Slj) 
for band 3. In particular, consistently with the analysis 
based on IC-NLS, band 2 has vanishing quasi- resonant 
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components ep2a,Spis, while the strongest components 
£pia, £p2s are amphfied thanks to the single decay pro- 
cess r2 + (— ri) — > (fl — (5w) + {—ft + Slu). In fact, the 
peak frequency Slo ~ 3.3 is given by the positive root of 
the equation Suj{Slu — 20) — 1 = 0, which is nothing but 
the nonhnear phase-matching condition of such process 
in normahzed form. Given the peak frequencies of such 
branches 2 and 3 {Suj ~ 3.3 and Suj ~ 6.2, respectively) 
the relative growing modes turn out to lie in the proxim- 
ity of the pump frequencies [as sketched in Fig. [TJc-d)], 
thus overlapping with band 1. Photons at one of such 
frequencies, however, unavoidably lead, according to the 
full NLS dynamics, to the generation of idler photons at 
mirror frequency, and both grow exponentially in pair 
(recall that they fall within band 1) at the expense of 
a single pump, as entailed by Manley-Rowe relations of 
the dominant process [see Fig. [TJb)]. For this reason we 
expect MI due to branch 1 to always prevail against the 
other two branches. 

This scenario has been validated by performing split- 
step simulations of Eq. ([T]) , using a two-pump input over 
a noisy background. The output spectrum at distance 
z = 20, displayed in Fig.[2][a) for il = 1.5, clearly shows 
the growth of noise over MI bands which are symmetri- 
cally located around all FWM orders. MI develops on top 
of the periodic exchange of energy between the pumps 
and FWM sidebands, as shown in Fig. [2jb). Maximum 
amplification occurs at a modulation frequency Slj ~ 0.9 
around all FWM orders, which agrees well with the peak 
gain frequency of band 1 [superimposed in Fig.[2ja)]. We 
find that the peak frequency 5uj is nearly constant with 
pump detuning f2, saturating to the value 6uj ~ I for 

> 3. At such larger pump detunings also the growth 
of harmonics of Su! become visible [Fig. [2][c)]. 

To further confirm that neither band 2 nor band 3 
leads to the exponential growth of asymmetric modes, we 
have repeated the simulation of Fig.[21[a) with additional 
seeds along the unstable eigenvectors of Fig. [TJc-d). The 
case of band 2 is displayed in Fig. [3lja): besides spon- 
taneous MI growing with the rate predicted by band 1, 
the output spectrum exhibits generated idlers with seed- 
idler pairs replicated around all FWM orders. As it is 
clear from in Fig. the idler is generated very rapidly 
from noise and then follows adiabatically the seed. Since 
the seed-idler pairs are then nearly symmetric and fall 
within Ml-band 1, they are observed to grow (on aver- 
age) with the gain characteristic of band 1 (i.e. with a 
gain slightly lower than the peak gain of spontaneous MI) 
and not following the much higher growth rate charac- 
teristic of band 2 (reported in Fig.[3l^b) for comparison), 
which proves the correctness of our interpretation. 

By extending our analysis also to the normal GVD 
regime, we have found that the low-frequency branch 1 
responsible for MI disappears. Our simulations of the 
NLS dynamics confirm that no MI occurs in this regime 
[see Fig. [2][d)]. Indeed the basic underlying mechanism 
of collective MI remains the decay of photon pairs from 
each pump into sideband pairs [see Fig. [Hb)], which 
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Fig. 3. (a) Output spectrum with additional input seed 
at Suj ~ 3.3 along the unstable eigenvector of band 2. (b) 
Power fraction evolution (in log scale) of seed (triangles), 
idlers (circles) , and peak of spontaneous MI growing from 
noise around the pump (thick solid line) and the n — 3 
FWM order (thin solid line). The dashed lines give the 
predicted growth rates associated with peak gain of band 
1 and 2 in Fig. (TJ respectively. 

can be nonlinearly phase-matched only in the anoma- 
lous GVD regime. Such modulation is then transferred 
via the Kerr effect onto the generated FWM orders. This 
explains why the modulation frequency Suj is fixed for all 
FWM orders and does not scale with the inverse square 
root of the relative power as would have been the case 
for a non-collective phenomenon ruled by conventional 
scalar MI. Moreover, here MI sidebands grow in an oscil- 
latory way [Fig.[3Jb)], owing to group- velocity mismatch 
between the pumps. 

In summary, we have found that FWM exhibits a col- 
lective MI process that leads to the exponential growth of 
sidebands around all multiple FWM orders, consistently 
with our Floquet analysis, once the latter is properly 
interpreted in terms of the associated eigenvectors. 
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